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The STAR collaboration at RHIC presents measurements of J/ip — > e"'"e~ at mid-rapidity and high 
transverse momentum {pr > 5 GeV/c) in p+p and central Cu+Cu collisions at a/s^ ~ 200 GeV. 
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The inclusive J /^jj production cross section for Cu+Cu collisions is found to be consistent at high pr 
with the binary collision-scaled cross section for p+p collisions, in contrast to previous measurements 
at lower pr, where a suppression of J/t/j production is observed relative to the expectation from 
binary scaling. Azimuthal correlations of J/ if) with charged hadrons in p+p collisions provide an 
estimate of the contribution of _B-meson decays to J/i/" production of 13% ± 5%. 

PACS numbers: 12.38.Mh, 14.40.Gx, 25.75.Dw, 25.75.Nq 



Suppression of the cc bound state J /ip meson produc- 
tion in relativistic heavy-ion coUisions arising from J /tp 
dissociation due to screening of the cc binding potential 
in the deconfined medium has been proposed as a signa- 
ture of Quark-Gluon Plasma (QGP) formation Mea- 
surements at ^ 17.3 GeV at the CERN-SPS ob- 
served a strong suppression of J/i/' production in heavy- 
ion coUisions [3] , although the magnitude of the suppres- 
sion decreases with increasing J/?/; px- This systematic 
dependence may be explained by initial state scattering 
(Cronin effect [31]), as well as the combined effects of finite 
J 1^ formation time and the finite space-time extent of 
the hot, dense volume where the dissociation can occur 

i- 

At higher beam energy (y^s^ = 200 GeV), the 
PHENIX collaboration at RHIC has measured J/ tp sup- 
pression for pt < 5 GeV/c in central (small impact pa- 
rameter) Au+Au and Cu-|-Cu collisions 0] that is similar 
in magnitude to that observed at the CERN-SPS. This 
similarity is surprising in light of the expectation that 
the energy density is significantly higher at larger colli- 
sion energy. It may be due to the counterbalancing of 
larger dissociation with recombination of unassociated c 
and c in the medium, which are more abundant at higher 
energy @,0,ii[3. 

Measurements of open heavy-flavor production may 
also shed light on J/ip suppression mechanisms. Non- 
photonic electrons from the semi-leptonic decay of heavy 
flavor mesons are found to be strongly suppressed in 
heavy- ion relative to p-|-p collisions at RHIC [HI, fl^ - 
an effect that has been attributed to partonic energy loss 
in dense matter [31 . This process may also contribute 
to high-pT J/ip suppression, if J/ip formation proceeds 
through a channel carrying color. 

The medium generated in RHIC heavy-ion collisions 
is thought to be strongly coupled [l3|, making accurate 
QCD calculations of quarkonium propagation difficult. 
The AdS/CFT duality for QCD-like theories may pro- 
vide insight into heavy fermion pair propagation in a 
strongly coupled liquid. One such calculation predicts 
that the dissociation temperature decreases with increas- 
ing J/ ip pt (or velocity) [15[ . The temperature achieved 
at RHIC (~ 1.5 Tc) [ij] is below this dissociation tem- 
perature at low J/tp Pt, and above it at > 5 GeV/c. 
Consequently, J/ip production is predicted to be more 
suppressed at high pT, in contrast to the standard sup- 
pression mechanism. This prediction can be tested with 
measurements of J/ip over a broad kinematic range, in 



both p+p and nuclear collisions. 

The interpretation of J/ip suppression observed at 
the SPS and by the PHENIX collaboration requires un- 
derstanding of the quarkonium production mechanism 
in hadronic collisions, which include direct production 
via gluon fusion and color-octet (CO) and color-singlet 
(CS) transitions, as described by Non-Relativistic Quan- 
tum ChromoDynamics (NRQCD) [31; parton fragmen- 
tation; and feeddown from higher charmonium states {xc, 
ip{2S)) and B meson decays. No model at present fully 
explains the J/ip systematics observed in elementary col- 
lisions [l7| . J/ip measurements at high-pT both in p+p 
and nuclear collisions may provide additional insights 
into the basic processes underlying quarkonium produc- 
tion. 

This letter reports new measurements by the STAR 
collaboration at RHIC of J/ tp production at high trans- 
verse momentum in p+p and Cu-t-Cu collisions at ^/s^ 
= 200 GeV [3|. The inclusive cross section and semi- 
inclusive J/i/j-hadron correlations are presented. 

The Cu-hCu data are from the RHIC 2005 run, while 
the p+p data are from 2005 and 2006. The online trigger, 
utilizing the STAR Barrel Electromagnetic Calorimeter 
(BEMC) [3| as well as other trigger detectors, required 
one BEMC tower with an energy deposition above a given 
threshold in coincidence with a minimum bias (MB) col- 
lision trigger The online trigger threshold, MB trig- 
ger condition, and sampled integrated luminosity for each 
dataset are listed in Tab. |T1 In Cu-|-Cu data, the most 
central 0-20% and 0-60% of the total hadronic cross sec- 
tion were selected as in [13, (111 ■ 

In this analysis, J/ip — > e+e^ (Branching Ratio 
(B)=5.9%) was reconstructed using the STAR Time Pro- 
jection Chamber (TPC) [13] and BEMC, with acceptance 
1 77 1 < 1 and full azimuthal coverage. Hadron rejection 
was achieved through the combination of BEMC shower 
energy, shower shape measured in the embedded Shower- 
Maximum Detector (SMD), and ionization loss (dE/dx) 



in the TPC [ll|,|23|. Electron purity is > 70% with high 
efficiency. At moderate px, the TPC alone can measure 
electrons with efficiency > 90% and sufficient hadron re- 
jection (- 10^) [Ulii]. 

Figure [1] shows di-electron invariant mass distributions 
for (a) p+p and (b) Cu-|-Cu collisions at y^s^ = 200 
GeV. The like-sign distribution measures random pair 
background from Dalitz decays and photon conversions. 
The J/ip mass window is 2.7 < Mf4 < 3.2 GeV/c^. 
Other correlated e~^e~ background is estimated to be < 
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TABLE I: Trigger conditions, off-line cuts and J/t/j signal 
statistics. Et is the BEMC trigger threshold. pTi and pt2 are 
the lower bounds for the two electron candidates. BBC (ZDC) 
means the coincidence of Beam Beam Counters (Zero Degree 
Calorimeters). S/B is the ratio of signal to background. 





p+p (2005) 


p+p (2006) 


Cu+Cu 


MB trigger 


BBC 


BBC 


ZDC 


Et (GeV) 


> 3.5 


> 5.4 


> 3.75 


Sampled int. lumi 


2.8 pb-^ 


11.3 pb-^ 


860 fib-'^ 


PTi (GeV/c) 


> 2.5 


> 4.0 


> 3.5 


PT2 (GeV/c) 


> 1.2 


> 1.2 


> 1.5 


J/iP PT (GeV/c) 


5-8 


5-14 


5-8 


J/ip counts 


32 ± 6 


51 ± 10 


23 ± 8 


S/B 


9:1 


2:1 


1:4 



(a) 200 GeV p+p 
5< Pt< 14GeV/c 




(b) 200 GeV Cu+Cu 
5 < Pt < 8 GeV/c 
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FIG. 1: (Color online.) Left: invariant dielectron mass distri- 
bution in (a) p+p and (b) Cu~f Cu collisions, for opposite sign 
(solid red) and same sign pairs (grey band) from data, and 
simulated J/ip peak for p+p (dashed). Right: J/tp pr distri- 
butions in p+p and Cu-|-Cu collisions at y^s^ = 200 GeV. 
Horizontal brackets show bin limits. Also shown are pertur- 
bative calculations for LO CS-I-CO (solid line) and NNLO* 
CS (band) direct yields, without feeddown contributions. 



10% 



H, [ill. Table U lists the offline cuts and J/ip 
signal statistics. Different thresholds were used for the 
two electron candidates, corresponding to different online 
trigger thresholds. 

The J/ip detection efficiency was calculated by two 
complementary methods. The first method was to de- 
termine the electron trigger efficiency by comparing trig- 
gered electron yield to the measured inclusive electron 
spectrum The non-triggered electron efficiency de- 
pends only on the TPC tracking efficiency, which was 
determined by embedding simulated electron tracks into 
real events [20|, and dE/dx efficiencies, determined from 
the distributions in real data [i^. The second method 
was to simulate J/ip events in PYTHIA [i^ , embed them 
into real events, and reconstruct the hybrid event to de- 
termine the J/ip trigger and detection efficiencies. The 



O ISR • STAR 
A UA2 o PHENIX 
V FNAL * UA1 
A CDF 



jrxlO^ 
(n=6.6) 




FIG. 2: XT distributions of pions and protons 1331. 1341 . 1351. 136 
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difference in estimated efficiency between the two meth- 
ods is < 10% for all datasets and is included into the 
systematic uncertainties of the inclusive spectra. This 
systematic uncertainty is correlated in p+p and Cu+Cu. 
A log- likelihood method is used to correct the J/ip effi- 
ciency and calculate the yields [29l |. 

Figure[T](c) shows the measured J/ip ^ e~^e~ pr spec- 
tra. The systematic uncertainties are dominated by kine- 
matic cuts, trigger efficiency (9%) and reconstruction ef- 
ficiency (8%), and are similar and correlated in p+p and 
Cu4-Cu. The normalization uncertainty for the inclu- 
sive non-singly diffractive p+p cross section is 14% [30| . 
Theoretical calculations shown in the figure are NRQCD 
from CO and CS transitions for direct J/ip's in p + p col- 
lisions [13 (sohd fine) and NNLO* CS result ^ (gray 
band). Neither calculation includes feeddown contribu- 
tions. The band for NNLO* gives the uncertainty due 
to scale parameters and the charm quark mass. The 
CS-I-CO calculation describes the data well and leaves 
little room for feeddown from ip' , Xc and B, estimated to 
be a factor of ^ 1.5. NNLO* CS predicts a steeper pT 
dependence. 

Proton and pion inclusive production cross sections 
in high energy p+p collisions have been found to fol- 

g{xT) / s^^"^ , where 



41|, |43: 



low xt scaling [4C 
XT ~ 2pt/\/s. In the parton model, n reflects the num- 
ber of constituents taking an active role in hadron pro- 
duction. Figure [2] shows the xt distributions of this data 
and previous J/ip, pion and proton data, from p+p coUi- 



271 . 138l . |39| cover the range 



sions. The J/ip data [25|, \2t 

^30 GeV to Vi =1.96 TeV. The J/ip exhibits xt 
scaling {n = 5.6 ± 0.2) at high pT, similar to the trend 
for pions and protons (n = 6.6 ± 0.1) [3, [s^]. While 
low PT J/ip production originates in a hard process due 
to the mass scale, subsequent soft processes could cause 
violation of xt scaling. At high pT, the power param- 
eter n = 5.6 ± 0.2 is closer to the predictions from CO 
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FIG. 3: (Color online). J/V' Raa vs. pr- STAR data points 
have statistical (bars) and systematic (caps) uncertainties. 
The box about unity on the left shows Raa normalization 
uncertainty, which is the quadrature sum of p-f p normaliza- 
tion and binary collision scaling uncertainties. The solid line 
and band show the average and uncertainty of the two 0-20% 
data points. The curves are model calculations described in 
the text. The uncertainty band of 10% for the dotted curve 
is not shown. 



and Color-Evaporation production (n ~ 6) [3l|, |43[ and 
much smaller than that from next-to-next-to leading or- 
der (NNLO*) CS production (n ~ 8) This is also 
evident from Fig. [T] (c). 

The nuclear modification factor Raa(j>t) 01 i defined 
as the ratio of the inclusive hadron yield in nuclear col- 
lisions to that in p+p collisions scaled by the underlying 
number of binary nucleon-nuclcon collisions, measures 
medium-induced effects on inclusive particle production. 
In the absence of such effects, Raa is unity for hard pro- 
cesses. 

Figure [3] shows Raa for JM vs pr, in 0-20% Cu-|-Cu 
colhsions from PHENIX and STAR, and 0-60% 
Cu-I-Cu from STAR. Cu-t-Cu and p+p data with pr > 
5 GeV/c are from STAR. The Raa systematic uncer- 
tainty takes into account the correlated efficiencies of the 
Cu-I-Cu and p+p datasets. Raa for J/i^ is seen to in- 
crease with increasing pr- The average of the two STAR 
0-20% data points at high-pr is Raa = 1.4±0.4 {stat.)± 
0.2 (syst.). Utilizing the STAR Cu-f Cu and p+p data re- 
ported here and PHENIX Cu-|-Cu data at high-pr SI 
gives Raa = 1.1 ± 0.3 (stat.) ± 0.2 (syst.) for pr > 5 
GeV/c. Both results are consistent with unity and differ 
by two standard deviations from a PHENIX measure- 
ment at lower pr {Raa = 0.52 ± 0.05 ^^). A notable 
conclusion from these data is that J/iJj is the only hadron 
measured in RHIC heavy-ion collisions that does not ex- 
hibit significant high px suppression. However, for the 
J/ip population reported here, the initial scattered par- 
tons have average momentum fraction ^ 0.1 (see also 
Fig. 2), where initial state effects such as anti-shadowing 
may lead to increasing Raa with increasing px- 

The dashed curve in Fig. [3] shows the prediction of 
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FIG. 4: (Color online). J/7/;-hadron azimuthal correlations. 
Lines show PYTHIA calculation of prompt (dashed) and B- 
meson (dot-dashed) feeddown contributions, and their sum 
(solid). 



an AdS/CFT-based calculation, in which the J/tp is em- 
bedded in a hydrodynamic model [i^ and the J/ip dis- 
sociation temperature decreases with increasing veloc- 
ity according to [ist . Its px dependence is at variance 
with that of the data. The dotted line shows the predic- 
tion of a two-component model including color screening, 
hadronic phase dissociation, statistical cc coalescence at 
the hadronization transition, J/^j formation time effects, 
and S-meson feeddown f^. This calculation describes 
the overall trend of the data. 

The other calculations in Fig. [3] provide a compari- 
son to open charm Raa- The solid line is based on the 
WHDG model for charm quark energy loss, with assumed 
medium gluon density dNg/dy ^ 254 for 0-20% Cu-fCu 



47| . The dash-dotted line shows a GLV model calcula- 



tion for D-mcson energy loss, with dNg/dy — 275 [4J 
Both models, which correctly describe heavy-flavor sup- 
pression in Au-|-Au collisions, predict charm meson sup- 
pression of a factor ^ 2 at pT > 5 GeV/c. This is in 
contrast to the J/^ Raa- This comparison suggests that 
high-pT J/'fp production does not proceed dominantly via 
a channel carrying color. However, other effects 0, |4^ 
may compensate for the predicted loss in this px range. 

Figure|4]shows the azimuthal correlation between high- 
Pt J/^' [pt > 5 GeV/c) and charged hadrons with 
PT > 0.5 GeV/c in 200 GeV p+p collisions. The J/^p 
mass window is narrowed to 2.9-3.2 GeV/c^ to increase 
the S/B ratio. There is no significant correlated yield in 
the near-side (A(f) ^ 0), in contrast to dihadron correla- 
tion measurements 15011. The lines show the result of a 
PYTHIA calculation [28[ , which exhibits a near-side cor- 
relation due dominantly to B ^ J/tjj + X. A fit to 
the data of the summed distribution (directly produced 
J/V' , feeddown from Xc, "0(25') and _B-meson) gives a 
contribution from _B-meson feeddown to inclusive J/^p 
production of 13% ± 5% at pr > 5 GeV/c. 

In summary, we report new measurements of J/^p pro- 
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duction in y/s =200 GeV p+p and Cu+Cu collisions 
at high pT [pT > 5 GeV/c) at RHIC. The J/ip inclu- 
sive cross section was found to obey xt scaling for px 
> 5 GeV/c, in contrast to lower pT J/tJj production. 
The J/ijj nuclear modification factor Raa in Cu+Cu in- 
creases from low to high pT and is consistent with no 
J/i/j suppression for px >5 GeV/c, in contrast to the 
prediction from a theoretical model of quarkonium disso- 
ciation in a strongly coupled hquid using an AdS/CFT 
approach. The two-component model with finite J/ip 
formation time describes the increasing trend of the J/tp 
Raa- Based on the measurement of azimuthal correla- 
tions and the comparison to model calculations, we es- 
timate the fraction of J/ip from _B-mcson decay to be 
13 ± 5% at PT > 5 GcV/c. 
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